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Abstract 

We study the conformal spectra of the critical square lattice Ising model on 
the Klein bottle and Mobius strip using Yang-Baxter techniques and the solution 
of functional equations. In particular, we obtain expressions for the finitized 
conformal partition functions in terms of finitized Virasoro characters. This 
demonstrates that Yang-Baxter techniques and functional equations can be used 
to study the conformal spectra of more general exactly solvable lattice models 
in these topologies. The results rely on certain properties of the eigenvalues 
which are confirmed numerically. 
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1 Introduction 

There has been much recent progress on understanding general conformal boundary 
conditions for rational theories on the cylinder [|l| and torus p| and their relation P] 
to integrable boundary conditions on the lattice. It is therefore of some interest to 
extend this understanding to other topologies such as the Klein bottle and Mobius 
strip. However, while much is known [|, ^, ^ about conformal partition functions 
in these topologies from the viewpoint of string theory, very little is known about 
integrable lattice boundary conditions on the Klein bottle and Mobius strip. Indeed, 
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Figure 1: Klein bottle (left) & Mobius strip (right). The boundaries are joined ac- 
cording to the orientations indicated by the arrows. The arrows at the bottom left 
corners on each face weight indicate its orientation. 



the lattice Yang-Baxter techniques which are well developed for the torus and cylinder 
have not yet been exploited in other topologies. 

In this paper, we use Yang-Baxter techniques to study the conformal spectra of 
the critical square lattice Ising model on the Klein bottle and Mobius strip shown in 
Figure 1. Although the Ising model has been studied previously in these topologies by 
free-fermion and Pfaffian p, |T0| techniques, we point out that these techniques 
do not generalize to other exactly solvable lattice models. In contrast, the methods 
based on Yang-Baxter techniques and the solution of functional equations developed 
in this paper will generalize to other exactly solvable lattice such as the A-D-E lattice 
models |]T3| . 

The layout of the paper is as follows. In Section 2 we review the Ising model 
on the torus. We define the representation of the Ising model and its single 
row transfer matrix T{u). We also summarize the solution [|l^, |l^ for the finitized 
conformal partition function using functional equations. In Section 3, we consider the 
Ising model on the Klein bottle. We define the double row transfer matrices D{u) = 
T{u)T{\ — u) and the flip operator F and discuss the simultaneous eigenvectors of F 
and D{u) and their relation to the eigenvectors of T{u). We calculate the conformal 
spectra of the Ising model on the Klein bottle in terms of finitized Virasoro characters 
by solving the relevant functional equation in the form of an inversion identity. In 
Section 4 we consider the Ising model on the Mobius strip. We define the double row 
transfer matrices D{u) with integrable boundary conditions corresponding to -|- and 
free boundaries. The relevant inversion identities for the conformal spectra of the Ising 
model on the Mobius strip are the same as for the cylinder with + or free boundary 
conditions on the left and right. For each eigenvalue, we identify the eigenvalue 
F = ±1 under the flip operator F according to the classifying pattern of 1- and 
2-strings in the complex w-plane. This identification is not obtained analytically but 
is confirmed by extensive numerics determining the parity F = ±1 of the associated 
eigenvectors. Finally, we end with a brief discussion. 
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2 Ising Model on the Torus 

In this section we briefly review the critical Ising model with periodic boundary con- 



ditions on the torus and its finitized partition function 12, O, M 



2.1 Ising model 

The Hamiltonian of the Ising model is given by 

^(W) = - E J^,cr^a,~hY,a, (2.1) 

<i,j> i 

where cxj = ±1 is the spin at site i, Jij > the interaction strength between site i and 
j, and h the external field. The sum is over all pairs of sites <i,j>. 

In this article, we consider only the nearest-neighbour critical Ising model on 
square lattice with zero external field. The partition function is given by 

Z{K,L) = ^exp(ir^ + a,cx,) (2.2) 

{o-} <i,j>h <i,j>v 

where K, L are the couplings in horizontal and vertical direction respectively. The 
model is critical on the line sinh(2i^) sinh(2L) = 1. The partition function is a sum 
over all possible configurations of spins on the lattice with/without constraints due 
to the geometry of the lattice and the boundary conditions. 

However, since the row-to-row transfer matrices in this representation do not com- 
mute, Baxter's commuting transfer matrices technique cannot be applied directly. 
Hence, we work in the ^3 representaion defined in the next sub-section for the follow- 
ing reasons: 

1. The critical Ising model on a square lattice can be conveniently formulated as 
the critical ^3 model. 

2. It can be easily generalised to other A — D — E models. 



2.2 representation and its transfer matrices 



The ^3 model is one of the Andrews-Baxter- Forrester models JT^ with L = 
introduced in 1984. In 1987, Pasquier constructed the critical A — D — E models |jl3 
which include the critical Al models. In the ^3 model, the spins a, b,c,d, . . . assigned 
to the sites of the lattice take heights from the set {1, 2, 3} and satisfy the adjacency 
condition that heights on adjacent sites must differ by ±1. 

Consider the Ising model on a square lattice of N columns and M rows where 
N=2L and M=2L' are both even. 
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Figure 2: Ising model couplings on A3 lattice. 



The Boltzmann face weights of the critical ^3 model are 
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where A = | is the crossing parameter and Sa = sin(aA) are the crossing factors. It 
is understood that the face weights vanish if the adjacency constraint a — 6 = ±1 is 
not satisfied along any edge. The spectral parameter u is usually taken in the range 
< u < A, so that the face weights are all positive but can be extended into the 
complex u-plane by analytic continuation. 

The ^3 model corresponds to the two dimensional Ising model. In ^3 model, the 
square lattice has two sublattices which can be even or odd. The heights a are even 
on the even sublattice and odd on the odd sublattice. On the even sublattice, the 
heights are fixed to the value 2. On the odd sublattice we identify the state a = 1 
with the usual + Ising state and a = 3 with the usual — Ising state. By the adjacency 
condition, if one of the sublattice is odd, then the other sublattice must be even and 
vice versa. The Ising model couplings act along the diagonals of the faces as depicted 
in fig. 0. It can be shown that the spectral parameter u is related to the couplings by 
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It is obvious that the interactions are isotropic at m = From this, we see that 
the ^3 model corresponds to the Ising model. Note that in the representation, 
it contains two mutually independent copies of Ising models. From now on, without 
stating explicitly, we use the term "Ising model" meaning "the Ising model in ^3 
representation" . 

The face weights ( |2.3| ) have the following local properties. They are symmetric 
under reflection about the diagonals 
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and satisfy the crossing symmetries 
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By ( |2.5| ) and ( p.6|) , we also have the following symmetries under horizontal and vertical 
reflections 
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The single row transfer matrix T[u) between the rows of heights a = {ai, 
and b = ■ ■ ■ , b^} has entries 



T(u 



a.b 
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(2.7) 

,aN} 

(2.8) 



with the periodicity a^+i = cti, ^tv+i = bi. 

Since the Boltzmann weights satisfy the Yang-Baxter equation and inversion rela- 



tion, the transfer matrices T{u) form a commuting family |jT^, |2T| with T{u)T{v) = 
T[v)T{u). Using ( |2.7| ) and the periodic boundary condition along the row we find 

T(A - n)„,5 = T{u\a = [T^('«)]a,5 (2.9) 

So T{X — u) = T^{u) and T{u) both belong to a commuting family of normal matrices 
which can be simultaneously diagonalized by a unitary matrix. The corresponding 
eigenvectors {Zn} are in general complex and independent of u. 

The adjacency matrix for the Ising model in the representation is 



A 




(2.10) 



i. e. Aij = 1 if heights i & j are admissible on adjacent sites and zero otherwise. It 
follows that the dimension of the transfer matrix is 

dim T(m) = TrA^ = 2(2^/2) _ 2^+1 (2.11) 

Note that the dimension is twice larger than the one in usual representation as defined 



in section p.l| since we are working in ^3 representation which contains two copies of 
Ising models. In fact, we can order all the row configurations into two blocks. In 
the first block the first height in each row is even. In the second block, we cyclically 
translate the row configurations of the first block by one unit to the right so that the 
first height in each row is now odd. Using the fact that T{u) commutes with the shift 
operator C = T(0) which is defined as 
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(2.12) 
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we obtain, in this basis, the block matrices 

where the orthogonal shift matrix C commutes with V{u) and each of the two blocks 
are of size 2^. It follows that the eigenvalues of T{u) occur in pairs ±T(n). As an 
aside, V{u) and CV{u) in ( p.l3| ), up to normalization, are equal to V{u) and W('u) 



as defined in chapter 7 of respectively. 



2.3 Finitized Ising partition function on the torus 

The critical properties of the Ising model are described by the unitary minimal con- 
formal field theory with central charge c = 1/2 and conformal weights A = 0, 1/2 and 
1/16. This conformal data is related to lattice quantities through finite-size correc- 



tions |T^, |T9[ to the eigenvalues of the transfer matrices. 

Consider the finite-size partition function Z^m of the critical Ising model on a 
periodic lattice of columns and M rows. The asymptotic behaviour of Znm in the 
limit of large N and M with the aspect ratio M/N fixed is given by 

ZT--(m) = Tr[Tiu)''] ~ exp [-NMf^Mu)] Z^^^'^^q) (2.14) 

where T{u) is the periodic row transfer matrix, /buik(w) is the bulk free energy and 
^•Torus^^-j is the universal conformal partition function with modular parameter 

g = exp(-27rie-^'"^). (2.15) 

Since the lattice is periodic, there is no boundary free energy. 

The finite-size corrections of the eigenvalues T{u) of the Ising model on a torus 
are given by []TB|, jTU 



iogr(M) = -Ar/b,ik + ^ 



f — - A - a) sin 4m - zfce"^™ + ifce^™ 
U2 / 



+ o(l) (2,16) 



where A, A are conformal weights and k and k are arbitrary non- negative integers 
yielding towers of eigenvalues above each primary level (A, A) with k = k = 0. The 
modular invariant partition function is expressed in terms of Virasoro characters as 

(g) = 5] Ar(A,A)g-^/i2+^+^^g'=^g^ (2.17) 

(A,A) k k 

= xo{q)xo{q) + xi/w{q)xi/i(i{q) + xi/2{q)xi/2{q) 

where q is the complex conjugate of q and the operator content is given by 

MA A) = / ^' = V2) (2 18) 

^ ' ^ 10, otherwise ^ ' ' 
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The c = 1 /2 Virasoro characters are given variously by 



oo 



= ^ n (i + V n (1 - ^'"^/') ^ (2-19) 

^ U=i k=i J 



1 

oo 



m>0 y-^rm j-_ 
771 even 



„-l/2 f oo oo ^ 

q-'"'\i/2{q) = '-l^\U{l + q'-'^')-U{l-q'-'^')\ (2.20) 



>fc=l k=l 



1 , 



m>0 
m odd 



(9)m (g)c 

g-'/'W6(g) = n(i + ^') (2-21) 



/c=l 

1 



;rm(»Tt+l) 1 oo 



m>0 



where {q)m = IYk=ii^ ~ 1^) for m > and (g)o = 1. The three different forms of each 
of the characters constitute the c = 1/2 Rogers-Ramanujan identities. 

For calculations, it is convenient to work with a finitized partition function nor- 
malized by the largest eigenvalue Tq{u) for which A = A = k = k = 







(1)" - 





Z^°™^(L;g) = V ^ |g|^/i2^^orus^^)^ AT, M ^ oo, M/A^ fixed (2.22) 



where n = 0, 1, . . . labels the eigenvalues Tn{u) of T{u) for N = 2L columns. The 
contribution |g|~'^/^^ from the largest eigenvalue Tq{u) in the limit M and large 
with M/N fixed is obtained in Appendix B by using the Euler-Maclaurin formula. 

Following |T^, |T6|, the eigenvalues T{u) of the Ising model on a torus can be 
obtained by solving the functional equation 

T{u)T{u + f ) = cos^^ 2u-R sin^-^ 2u (2.23) 



where L = N/2, T{u) = T(f - u) and T(0) is the shift operator. Here = ±1 is the 
eigenvalue of the spin or height reversal operator 

N 

Ra,b = l[KA-b, (2-24) 
which satisfies = I, RT{u) = T{u)R. For R = 1 the solution of the functional 



equations is ||T6| 

T(n) =ev^(2e2^"+-'/^)""n (e'" + ^/^fctan((^)) (2.25) 



-L ^ 



k=l 
7 



with = /i^ = 1 for all k and 11^=1 /^fc = 1- Similarly for i? = — 1, the solution is 

L-l 
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(2.26) 



with = yu| = 1 for all k. Up to an overall sign, the eigenvalues are thus determined 
by the set {^k}- We note that the total number of eigenvalues is 2-^ + 2-^ = 2-^+^ and 
that each eigenvalue appears with its negative. 

It is often desirable to remove this trivial degeneracy of the eigenvalues ±T{u). 



Following |T6[ it is convenient to choose 
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e= n /^fc' 
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e = n l^k, R 

fc=L{L+2)/2j 



(2.27) 



This leads to the finitized Ising torus partition function [ITB]: 
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(L;g) = Xo(LtJ;g)Xo(L^J;g) + |g|Xi/2(LtJ;g)Xi/2(L^J;g) 

+ |gr/'^i/i6(L^J;g)Xi/i6(LfJ;?) 



with finitized Virasoro characters 120 
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(2.30) 



n (i + ^'-^/') - n (i - ^ 



k-l/2 
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m — 1 
m odd 



(m2-l)/2 
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E n 1-^(1 
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]^(l + g^) = g™(™+l)/2 
fc=l m>0 



(2.31) 



Here {yuji denotes the sequences /ii,...,/^/, with /i^ = ±1 and {/i}^ denotes the 
subset of sequences with 11^=1 /^fc = =tl- The Gaussian polynomials or g-binomials are 
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m. 



L—m Ii 



Ml 



E E---E^^^ 

/l=0/2=0 lm=0 



h + .-.+Ir, 



{g)rn{q)L-m. ' 

0, 



<m< L 
otherwise 



(2.32) 



with (g)m = n™=i(l ~ The R = 1 sector of the spectra is given by Xq and X1/2 
and the R = —1 sector is given by Xi/iq. 
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3 Ising Model on the Klein Bottle 



3.1 Flip operator and partition function 

To form a Klein bottle we consider an x M lattice of spins j and impose peri- 
odic boundary conditions along the horizontal direction of the lattice by identifying 
column + 1 with column 1 

ai,Ar+i = ai,i, i = l,2, ...,M. (3.1) 

Along the vertical direction, we identify row M + 1 with the left-right flip of row 1 

«M+l,j = 0,l,N+2-j, J = 1; 2, . . . , A^. (3.2) 

This flip reflects about the column j = 1, or equivalently j = ^ + 1, so that the spins 
at j = 1 and j = y + l = L + l are left invariant and the even and odd sublattices 
do not mix under the flip. 

The flip operator F that implements the flip ( p.2|) has entries 

N 

Fa,b=Y[ ^a,,bN+2-j (3-3) 

This matrix has the same dimension as the transfer matrix T(u) and satisfies -F^ = I, 
F'^ = F so that F has eigenvalues F = ±1. Clearly the flip operator F can be 
expressed, in a suitable basis, as a simple block matrix 



(3.4) 



where In is the n x n identity matrix, r is the number of row configurations which 
are invariant under the flip F, 2s is the number of non-invariant row configurations 
which are related in pairs by the flip F and r + 2s = dimF = 2(2^/^) ^ 2^+^ 

The number of left-right symmetric or invariant row configurations is simply the 
number of configurations of length N/ 2 and so 

r=y\A^/^] , 2s = TrA^- y ^3^5) 

^ L J a,fe ' ^ L J a,b ^ ^ 

a, 0=1 a, 0=1 

Diagonalizing the matrix A and using spectral decomposition we find that 

^ / 1 \' N / 1 \' f 3(2^/2), L = f even 

V ' ' V I 2(^+3)/2^ L = f odd 




9 



We see immediately that F has r + s eigenvalues +1 and s eigenvalues —1 and so r 
is the net number of positive eigenvalues 



(3.7) 



The flip operator F does not commute with the single row transfer matrix T{u) 
so it is not possible to simultaneously diagonalize T{u) and F. Instead, we find 
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[Ti\-u)F], 



T^(u)F 



(3.8) 



a,b 



For this reason we introduce the double row transfer matrix D{u) defined by 

D{u) = T(A - u)T{u) = T^{u)T{u) (3.9) 

Clearly D{u) is real symmetric and positive definite for real u. Moreover, F and 
D{u) do commute so they can be simultaneously diagonalized 



FD{u) = FT{\ - u)T{u) = T{u)FT{u) 

= T{u)T{X-u)F = T{X-u)T{u)F = D{u)F 

The partition function of the Ising model on the Klein bottle is 



7Klein , 
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(3.10) 



(3.11) 



where n labels the eigenvalues with respect to the common set of eigenvectors of 
F and D{u). Although the fiip operator F refiects about the column j = 1 we could 
have chosen a fiip operator Fk that refiects about the column j = k. In this case 



Fk = C^FC 



(3.12) 



It is then easy to check that the partition function (|3.11|) is translation invariant by 
using the cyclicity of the trace and the fact that C and D{u) commute 



Tr 



FkDiu 



.M/2 



Tr 
Tr 



C^FC-^D{u)^^'^ 
FDiut"^ 



Tr 



FC-''D{uY'I^C^ 



(3.13) 



Hence we only need to work with the fiip operator F = Fi defined in ( p.3| 
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We will show that the simultaneous eigenvectors of F and D{u) are 



'1 5 



-1 ! 



(3.14) 



with 



Fx, = +x„ j = l,2,...,r; Fz^ = ±zt, k = l,2,...,s (3.15) 

so that FD{uY^^'^z^ = ±D{u)'^^/'^zf. It follows that the eigenvalues corresponding 
to z'l and 2^ cancel in ( p.ll|) . So only the eigenvalues of D{u) corresponding to 
the even eigenvectors {xi, • • • , Xr} contribute to the partition function. We will show 
that the eigenvectors {xj} are precisely the real eigenvectors of T{u) whereas the 
eigenvectors {2^} are simply related to the complex eigenvectors of T{u). 



3.2 Simultaneous eigenvectors of D{u) and F 

In this section we study how the simultaneous eigenvectors of D{u) and F are related 
to the eigenvectors of T{u). Consider the transfer matrix T{u), with u in the closed 
interval [0, A] C R. This is a commuting family of real normal matrices that along 
with the flip operator F satisfies 

T{u)T{v) =T{v)T{u), T^{u) =T{X-u), FT{X - u) = T{u)F 

D{u) = T^{u)T{u), FD{u) = D{u)F F^ = I (3.16) 

The entries of T{u) are real analytic functions of u that can be analytically contin- 
ued into the complex plane u E C. In general, the eigenvalues T{u) and eigenvectors 
{Zn} of T{u) are complex with the eigenvectors z G {Zn} being independent of u. 
Specifically, from ( [A.l| ) we have 

T{u)z = T{u)z, T^{u)z = T(A - u)z = T(A - u)z = T{u)z (3.17) 

The expressions for T{u) can be obtained in the following way. Given an eigen- 
vector z = {z^, z^, . . .) with constant complex entries (since z is independent of u), 
we expand the product T{u)z along the row with the largest absolute value entry 
in z. Since z is an eigenvector of T{u), the corresponding eigenvalue T{u) is 

^(«) = 4E[T(«)k,^^' (3.18) 
^ j 

Since T{u) is a finite matrix with entries being analytic functions, it follows that all 
the eigenvalues T{u) which, are the finite sums of analytic functions, are analytic 
functions of u. 

We say that two eigenvalues are degenerate if they agree on the whole interval 
[0, A]. Clearly, by analytic continuation, two eigenvalues are degenerate if they agree 
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on any open interval in [0, A]. Consequently, non-degenerate eigenvalues can only 

cross at isolated points. The difference between these non-degenerate eigenvalues is 
an analytic function with isolated zeros which cannot accumulate in the closed interval 
[0, A]. It follows that there is only a finite number of points where the crossings of 
the various eigenvalues can occur. Therefore at a "generic" point u — Uq there are 
no accidental crossings, that is, any two eigenvalues are either degenerate on [0, A] 
or distinct at -u = uq. We conclude that the simultaneous eigenvectors of T{u) are 
uniquely determined (up to linear combinations for the degenerate eigenvalues) by 
the distinct eigenvectors at u = uq. So, in practice numerically, it is sufficient to fix u 
to a "generic" value Uq and diagonalize T{uo) to find the simultaneous eigenvectors. 

Let {Zn} be a set of orthonomal complex eignvectors of T{u). We want to construct 
a set of simultaneous eigenvectors of D and F. Suppose, after 

suitable normalization, that the eigenvector z = x of T{u) is real. Then it follows 
that the corresponding eigenvalue T{u) is real for all real u. We refer to these as the 
real eigenvalues of T(u). 

T{u)xj = T^{u)xj = T{u)xj, D{u)xj = T{ufxj, j = 1, 2, . . . , r (3.19) 

In the remaining case we have to work with complex eigenvectors with both real and 

imaginary parts non-zero. In this case, the eigenvalues appear in complex conjugate 
pairs since z and z are linearly independent and T{u)z = T{u)z implies 

T{u)z = T{u)z = T{u)z = T{u)z. (3.20) 

Hence 

T{u)zk ^ T{u)zk, T'^{u)zk^T{u)Zk, D{u)zk ^ \T{u)\''zk, A;=l,2,...,s 

(3.21) 

with a similar statement for z^. 

Let {xi, . . . ,XriZi, . . . ,Zs\ be the set of mutually orthonormal real and complex 
eigenvectors corresponding to the real and complex eigenvalues T{uq) G R, T{uq) G C 
with ImT(Mo) > 0. Then the r + 2s mutually orthonormal simultaneous eigenvectors 
of D{u) = T^{u)T{u) and F are 

{a?!, . . . , Xr, x'^ , . . . , a;^, x^^ , . . . ,Xg} (3.22) 

where 

Fx^^Fx^, xt^zt, x^^-izl z^ = ^{I±F)zk (3.23) 

and the corresponding eigenvalues of {T(A — u)T{u), F} are {T{u)\ F}, {\T{u)\^, +1} 
and {|T(m)P,— 1} respectively for the three groups. Although it is not apparent 
here we will see later that, since Fzk — Zk, the vectors Xj, a?^, x'j^ are all real with 
x'^ = V2Re(zk), ajfe = V2Im(zk). 
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We show that ( |3.22| ) is an orthonormal set of eigenvectors and hence a complete 
set of hnearly independent eigenvectors of T(A — u)T{u) and F. To estabhsh this we 
note some properties of T{u) and F (see Appendix A): 

• Two eigenvectors 2j and Zj of T{u) are orthogonal if their eigenvalues are dis- 
tinct at M = Uq. 

• If T{u)z = T{u)z then Fz is also an eigenvector of T{u) with the corresponding 
eigenvalue T{u) since 



T{u) (Fx) = {T{u)F)x = (FT(A - u)) x = T{u)Fx (3.24) 

Now we observe that the eigenvalues corresponding to the eigenvectors Zk and Fzi are 
distinct since they have imaginary parts ImT(n) > and lm.T{u) < respectively. 
So these two eigenvectors are orthogonal 

zlFzi = 0, k,i = l,2,...,s (3.25) 

The required orthonormality now follows by straigthforward calculation of inner prod- 
ucts 

(4)^ = -^ziiI±F)x, = ^il±F)zlx, = (3.26) 
{z^^zf = \zi{I±Ffz, = zl{I±F)z, = zlz, = 5,,, (3.27) 

= \zI{I±F){ItF)z, = Q (3.28) 
Similarly, the required eigenvalues follow by straightforward calculation 



T(A - u)T{u)xj = T{u)T{u)xj = T{uYxj, Fxj = Fxj (3.29) 
T(A-n)T(«)4 = -^TiX-u)Tiu)iI±F)z, 

1 



^ T(A - u) [T{u)zk ± T{u)Fzk) (3.30) 



^ {\nu)\'z,±\T{u)\'Fz,) = \Tiu)\ 



2„± 



V2 



Fzt = -^F{I±F)z, = ±^{I±F)z, = ±zt (3.31) 
3.3 Eigenvectors of T{u) 

To study the simultaneous eigenvectors of T(u) and T'^{u), we fix m = uq to a 
generic value to avoid accidental crossings of eigenvalues. We also set T = T{uq) 
and T = T{uq) so that the eigenvector equations for the common eigenvectors are 

Tz = Tz, T^z = Tz, Dz = T^Tz = \T\'^z (3.32) 
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The transfer matrix T does not commute with the flip operator F but the double 
row transfer matrix D = T^T does. In this section we obtain the structure of the 
complex eigenvectors z. 

Let us order the r + 2s paths a into three blocks of size r, s and s respectively. In 
the first block we put the r paths that are invariant under the flip F. The remaining 
2s paths occur in pairs related by the flip and we place these in order in the second 
and third blocks. In this basis, the block structure of F, T and is 

/ A 

G and we have used the fact that FT 
The r + s even and s odd real eigenvectors of F are of the form 





(I 









(A B 


F = 








;) 


T = 


G^ D 




lo 


/ 






[b^ G 


where A'^ = 


A, 




= H andG'^ = G and 




(3.33) 



T^F. 





( ° 




c 







(3.34) 



It follows that the r + s even and s odd common eigenvectors oi D = T'^T and F 
are real linear combinations of these r + s even and s odd eigenvectors respectively. 
Since the real symmetric matrices D = T'^T and F commute it follows that the real 
and imaginary parts of z in (|3.32|) are common real eigenvectors of D and F of the 
form( pIMD . Consequently, there are three possibilities (i) z is even and Fz = z, (ii) 
z is odd and Fz = —z or (iii) z is mixed so that after normalization Fz = z (In 
fact, it is not hard to show that Fz = e^^z for any arbitrary real 6' if z is mixed). But 



Tz = Tz 



T^z 



T^iFz 



(T 







T{Fz) 
Im(r) z\Fz) 







(3.35) 



(3.36) 



Tz and 
-T)z\Fz) = 

So if T ^ R, we conclude that z and Fz are orthogonal. Hence Fz ^ ±z and 
consequently Fz = z and the complex eigenvectors of T occur in complex conjugate 
pairs 

a \ f ^ 

b + ic \ , Fz = z = b — ic 
b — ic J \b + ic^ 

We can write the eigenvector equations for T and T'^ in terms of the symmetric 
and skew-symmetric combinations 

B + G 
D + D^ 
H + G 

B-G 
D-D'' 
G-H 



(T + T')z 



{T-T')z 





2Re(T)z (3.37) 



2zIm(T)z (3.38) 
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We now look for a complex eigenvector of the form ( ^.36 ) and take the real and 
imaginary parts in ( p.37] , p.38|) 




2 Re(T) 



/ 

2 Re(T) c 



-2Im(T) 



2 Im(T) 



a 



(3.39) 



(3.40) 



(3.41) 



(3.42) 



The symmetric equations can be simplified and from the two skew-symmetric equa- 
tions we obtain two eigenvector equations for real symmetric matrices 



2A 

V2{B + Cf 



V2iB + C) \f a \ 
D + D^ + H + GJ \V2bJ 



H 



-G 



H + G 
D + D^ 



2 Re(T) 




-4(Im Tf 



-4(Im Tf 




(3.43) 
(3.44) 

(3.45) 
(3.46) 



Equation ( |3.4(j| ) now reduces to a simple eigenvector equation for a real symmetric 
matrix 

(3.47) 



[2{B - GY {B-G)-{D- D' Y + (H - GY]c = 4(Im T)'c 



From numerical observation, this equations yields s real eigenvectors with Im T ^ 
at any generic point u. This is verified in section p.5| by counting the number of 
real eigenvalues of T which matches the net number of positive eigenvalues of F in 
(p.6|). Indeed, we will see that this generic situation leads to s complex conjugate 
pairs of eigenvectors of the form ( p.36| ). A non-generic choice of u would lead to fewer 
than s complex conjugate eigenvalues and eigenvectors. From ( |3.44| ) we see that. 
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in the generic situation, the s vectors c must also satisfy the simple real symmetric 
eigenvector equation 

[D + - H -G]c = 2 Re(r)c (3.48) 

For a given solution vector c we have thus determined the real and imaginary parts of 
the eigenvalues T = Re(T) ± zlm(T). The vectors a and b are determined by ( |3.42D 

a = (Im T)-\B-C)c, 6 = (2 Im T)-\D-D^ -H+G)c, lm{T) ^ (3.49) 

We have thus completely determined the s pairs of complex eigenvectors z, z of ( p.36| ) 
corresponding to the complex eigenvalues T = Re(T) ± ilm{T) with Im(T) ^ 0. 
The remaining equations such as ( |3.43 ) that we did not use in deriving the these 



eigenvectors must be automatically satisfied. 

Each of the s complex conjugate pairs of eigenvectors gives rise to one even eigen- 
vector 33+ = v^R-e(z) and one odd eigenvector x = \/2lm(z) of F and D = T'^T. 
These completely exhaust the odd common eigenvectors of F and D = T'^T. In the 
generic situation, it therefore remains to obtain the remaining r real even eigenvectors 
of T 

fa\ 

x=\b\ (3.50) 

which necessarily satisfy 

Fx = x, Tx = T'^x = Tx, Dx = T'^Tx = \T\'^x, T G R (3.51) 

These eigenvectors are in fact given by the symmetric eigenvector equation ( p.43|) 
which yields r + s even eigenvectors. The r solutions that we need are the solutions 
that remain after the s solutions for (a, b) obtained in the complex case are removed. 

In a non-generic situation, such as u = when T reduces to the shift operator, 
there are fewer than s complex conjugate pairs of eigenvectors of T. In this case, 
each complex conjugate pair that is lost is replaced by one even and one odd real 
eigenvector. This is the only time that odd real eigenvectors of T occur. We therefore 
conclude that, in the generic situation, all eigenvectors z of T satisfy 

Fz = z (3.52) 

In this case the common real eigenvectors of D = T'^T and F are 



X 



X 



1 

-^{I + F)z = V2\b\= V2Re{z) (3.53) 
V2 [b) 



-^{I ~ F)z = V2 c I = v^Im(2) (3.54) 



V2 

Clearly, these eigenvectors are of the form ( p. 34 ). 
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3.4 Conformal finite-size corrections 



Consider the finite-size partition function Zjmm of tlie critical Ising model on a Klein 
bottle formed from a lattice with columns and M rows. The asymptotic behaviour 
of Zi<^M in the limit of large N and M with the aspect ratio M/N fixed is given by 

Zf^^u) = TT[FD{uf'/'] ~ exp [-NM^^,M ' M/bdy(w)] Z'^'^'^iq) (3.55) 

where D{u) is the double row transfer matrix, /buik(w) is the bulk free energy, /bdy('w) 
is the boundary free energy and Z^'^™(g) is the universal conformal partition function. 
Since we are working with a double row such that D{u) = T{u)T[\ — u) and that 
9Torus(A — u) = qroTusiu), the product qTomsiu) gTorus(M) is real positive and thus we 
define the modular parameter to be real and given by 

q = exp(^- 2-71 — smAuj = |gTorus| (3.56) 

where gxorus is given by ( p.l5| ). Because of the presence of the flip, the left and right 
movers mix. Consequently, as is the case on the cylinder, there is only one copy of 
the Virasoro algebra. The leading corrections to each transfer matrix eigenvalue D{u) 
are obtained from the torus case by identifying A with A and k with k in (|2.16| ) 



1 Air / c \ /1\ 

- log D{u) = -iV/buik(u) - /bdy(u) + ^ ( ^ - ^ - ^) sm(4n) + o (-) (3.57) 



with A = 0, 1/2, 1/16 and k a non-negative integer. Cardy argued that as the 



model is conformal invariant, all spectra are organized into Virasoro characters, and 
thus the eigenvalues naturally divide into three towers with 

^Klein^^) = EA^(A) (g^)-'=/^'+^E^'' = EA/'(A)XA(g^) (3.58) 
A fc A 

where A/'(A) G {0, 1,2,.. .} is the operator content on the Klein bottle. Equation 
( p.58| ) with Af{A) = 1 is confirmed numerically in the next section. 



3.5 Partition function of the Ising model on the Klein Bottle 

The finitized conformal partition function on the Klein bottle is 

Z'''''^L;q) = j:Fni^) (3.59) 



where n = 0, 1, . . . labels the eigenvalues D„ of D{u) and F„ of F for N = 2L 
columns and Dq is the largest eigenvalue. We will obtain this partition function from 
the eigenvalues T{u) of the periodic row transfer matrix as given by ( |2.25| ), ( |2.26| ). To 
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do this we need to separate out the r real eigenvalues from the 2s eigenvalues which 
appear in complex conjugate pairs. 



It can be shown |^ that the eigenvalues T(u) possess the periodicity 



T(« + -) = (-!)- TH 



(3.60) 



and that within the period strip — | < Re{u) < in the complex w-plane the 
eigenvalues have exactly 2P zeros where 

L, R = +l 
R=-l 



P 



(3.61) 



For the ground state, given by /x^ = 1 for all k, the zeros occur within the period strip 
as 2-strings on the lines Re(M) = ~f 5 For excitations, zeros also occur on the lines 
Re(M) = l,-^- The eigenvalues are in fact completely determined by the patterns 
of zeros in the complex u-plane. Indeed we will show that T{u) can be uniquely 
expressed in the form 



T{u) = Y]_ sin2(M — m^), R^, mi, ■ ■ ■ , wp G C 



(3.62) 



k=l 



and obtain explicit expressions for the constant i?^ and zeros Uk- 
Let us define the monotonically increasing sequences 



tan 



tk 



tan(ff 



R 
R 



so that 



tk — tp+i-kJ ^ 



I 2 



P odd; 



k= 1,...,L 

k = l L 



1,...,P 

p 

k=l 



(3.63) 

(3.64) 
(3.65) 



For R = +1, we can now re-express (|2.25| ) in the form 
T{u) = ev^(2e^»+x)-"n(e4- + ,;,,tan((^) 



k=l 



(3.66) 



k=i 



g2i(«+|logtfe-fMfe + f ) _ g-2i(«+ilogtfc-f/Xfe + f ) 

eV2e'?[-^+^^=i'^^] n sin2(M + i logt^ - f/i^ + f 

k=l 

P 

R+ W sin2(M - Uk) 
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k=i 



where = /i^ = 11^=1 /^fc = 1 and we have used ( p.65| ). Suppose, for a given eigenvalue 
T{u), that m is the number of fXk such that fXk = —1- Since 11^=1 /i/c = 1, '"^ = 2p is 
even and 



k=l 



fik = + {L — m) = L ~ 2m = L — Ap 

Hence 

so that ( p.66|) becomes 

T(u) = 
where 



sin2(M — 



1(5 



L mod 2,1 



LL/2J 



(3.67) 
(3.68) 



nsin2(M-M+)sin2(M-n^) (3.69) 

A;=l 



TT TT 2 

-^/ifc - T - 7 log^fc, 



4 4 



A; = 1,2,...,P 



Ml 



Mi 



TT TT 

Mfe = -/ifc - ^ - ^logtfc, A; = 1,2, [P/2J 

TT TT 

= -pk - - + - logtfc, /c = 1, 2, . . . , [P/2J 



and 



(3.70) 
(3.71) 
(3.72) 

/ifc = /ip+i-fc, /c = 1, . . . , [P/2J. (3.73) 

Here u^T are the zeros in the upper and lower half planes respectively. Note that, 
if pk = fik, then is simply the complex conjugate of u^. Note also that real zeros 
Mfc e R occur for 

A;=[P/21, Podd (3.74) 
Similarly, for R= —1, ( [^.26| ) becomes 

L-l 



T(u) 



n (e^'^^ + ^/Xfctandf 



k=l 



L-l 



fc=i 



(3.75) 



L-l 



fe=l 



e v/Le^(^-^)e'?S^=i ^'^ H sin2(M + f logt^ - f/i^ + f ) 
p 

P~ Y[ sin 2(-u — Uk) 



k=l 



where = yu| = 1 for all k and we have again used (|3.65|) . But now 



exp 



T 



L-l 
k=l 



exp 



( 2m 

V A 



(3.76) 
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where m is the number of fik such that jj^ = 



— 1. Hence ( ^.751 ) reduces to 



T(u) =e(-z)™yL[sin2(M-MpL^.)j^^"°'''° H sin2(M - <) sin2(u - ) (3.77) 

' k=i 

We want to pick out the eigenvalues of T{u) which are real for all m G R. First, 
suppose the zeros occur in complex conjugate pairs and is real. Then it is clear 
that T{u) is real for all u G R. Conversely, suppose T{u) is real for all m G R. Then 

p 



T{u) =R^Y[sm2{u-Uk) = T{u) , mgR (3.78) 



k=l 



which implies that the zeros Uk occur in complex conjugate pairs and that is real. 
Thus, T(m) is real for all m G R if and only if R^ is real and the zeros Uk occur in 
complex conjugate pairs, that is, if the patterns of zeros in the upper and lower half 
planes are related by complex conjugation. Since A and k are determined by the 
pattern of zeros in the upper half plane and A and k are determined by the pattern of 
zeros in the lower half plane, this picks out precisely the eigenvalues for which A = A 
and k = k |TB|. By ( |3.69| ) and ( p. 771 ), we see that the eigenvalue T{u) is real for all 



M G R if and only if 



fik = /ifc and Y[fik = l. (3.79) 

k=l 



For R = —1 and L even (P odd), the last condition in (|3.79|) reduces to fi^L^-^ = 1. 

Using ( |3.79| ) and allowing for e = ±1, we can count the number r of real eigenvalues 
T{u) G R 

' 3(2^/2N I 

) 3.80 

2(L+3)/2 ^Q^^ ^ > 



2(2L#J +2LVJ^ = I 



in agreement with (|3.6|) . Since r is the net number of positive eigenvalues of the flip 
operator F, as discussed in Section 3.1, any real eigenvector of T{u) is an eigenvector 
of F with F = +1. The remaining 2s eigenvectors occur in complex conjugate pairs 
with F = ±1. After cancellation of the contributions from the complex eigenvalues 
and removal of the trivial ±T(u) degeneracy, the partition function for the Ising model 
on the Klein bottle is 

«"(«) = T.FnD4u)''/' = ^'jz Dn{uY"^ = 'f:' \T^{u)f (3.81) 

n n=0 n=0 

where the sum is over the r/2 real eigenvalues of F){u) 

Dn{u) = \Tn{u)\\ n = 0,l,...,r-l, T„(m) G R (3.82) 
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3.6 Finitized Ising partition function on Klein bottle 

In this section, we follow [|1^ to obtain the finitized conformal partition function 
^-Kicinj^^. terms of finitized Virasoro characters. 

We begin by using ( p.27| ) to remove the degeneracy in the eigenvalue expressions 
(p.25|) and (|2.26|) for the single row transfer matrix eigenvalues T{u). This degeneracy 
is irrelevant for the eigenvalues D[u) = T{u)T{u) of the double row transfer matrix. 
Using ( ^1271) and (|2:25| ) we find for = +1 



D{u) = T{u)T{u) 

LL/2J 



-2L+1 



e + i/ipL-i tan( 



- 7r(2rL/2l-l) - 
4L 



2<5r 



X 



n e^^" + z/ifctan(^i^) /i,e^^" + ^cot 



k=l 



'7^(2fe-l)^ 
^ 4L ' 



(3.83) 



with Hk = HL-k+i- Similarly, from (|2.26|) we find for R 



D{u) = T{u)T{u) = 2^(1-^)^ 

X n 



,4m 



tan( 



2L 



25 



L mod 2,0 



e^'" + ^/ifetan(ff) ^x^e^- + ^ cot(ff ) 



k=l 



(3.84) 



with /ifc = /iL_fc. We remark that although the forms for T{u) are not unique due to 
the choice of e, the forms given here for D[u) are unique since they are independent 
of e. 

Taking the ratio of ( |3.83| ) and (|3.84| ) with the largest eigenvalue Dq{u) given by 
( |3.83D with all fik = fik = ^ and using the limits 



lim A^log 



itan(5f 



o4m 



+ i tan 



lim A^log 



_^iiu _|_ ^ i^Q^ 



' ttK •< 

, AT y 

• ttK ' 

, AT . 



o4i?i 



+ ^cot(5f) 



-2mire-^*'' 
27riire^^'' 



(3.85) 



along with the condition ( |3.79| ) and the definition of q in ( p.56| ), we obtain 
Z^'<='°(L;g) = 

L(L-1)/2J 



LL/2J [1/2] 



+igr/ 



E n ^'''-'-^ 



(3.86) 



Here the factor \q\^^^ arises from the scaling limit of the ratio of the largest eigenvalue 
in ( p.84|) , which has all /Xfc = = 1, to the overall largest eigenvalue Dq{u) in ( p.83|) 
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which has all Hk = fik = 1- This ratio can be computed using the Euler-Maclaurin 
formula as in Appendix B. Using the expressions for the finitized Virasoro characters, 
we finally obtain the finitized partition function for the Ising model on the Klein bottle 

Z^'--(L;g) = Xo(LfJ;g') + |g|Xi/2(LfJ;g') + |gr/«Xi/i6(L^J;g^) (3.87) 

Notice that if we set q = I, then Z{L; 1) = r/2 as given by (|3.80| ) corresponding to 
the counting of the number of real eigenvalues left after cancellation of the complex 
conjugate pairs. As for the torus, the even sector of the spectra with i? = 1 is given 
by Xq and X1/2 and the odd sector with i? = — 1 is given by Xi/iq. 



4 Ising Model on the Mobius Strip 

In this section we consider the Ising model on the Mobius strip. Consider a square 
lattice consisting of a strip of N columns and M rows. Along the vertical direction we 
apply a flip operator F before joining row M to row 1. We apply integrable boundary 
conditions E^ on the left and right edges of the strip. These boundary conditions 
labelled by +, — , F (free) are specified by boundary triangle Boltzmann weights on 
the left and right edges of the double row transfer matrix D{u). We build up the 
Mobius strip by M applications of the double row transfer matrix. For compatibility 
with the flip the boundary conditions on the left and right must be the same. This 
is because the Mobius strip only has one edge. By the spin reversal symmetry, the 
partition function with + boundary conditions is the same as the partition function 
with — boundary conditions so we do not need to consider — boundary conditions. 
The left and right boundary weights are 



B 



u 




B b 



u 




(4.1) 



where it is understood that these weights vanish if the adjacency condition is not 
satisfied along the edges (a, h). These weights must satisfy the boundary Yang-Baxter 
equation. For the Ising model, we study the boundary conditions specified by the 
following non-zero weights 



Fixed + : The edge spins are fixed to a= 1 

•1 



B 



u 



B 2 



u 



(4.2) 



Free F : The edge spins are fixed to a=2 



B 



u 



Be 



u 



1 

71' 



1,3. 



(4.3) 
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For a strip with columns, the double row transfer matrix D{u) with rows 
a = {ao, . . . , Oat} and b = {bo, . . . , b^} is defined diagrammatically ||2l| by 




For the Ising model on the Mobius strip we require aQ = b^ = = b^j. The double 
row transfer matrices are real symmetric and crossing symmetric 

D{u) = D'^{u) = D{X-u), MGR. (4.5) 

The dimension of the transfer matrices D{u) are given in terms of the adjacency 
matrix A by 

ATI f2^-\ fixed + 



dim£)(M) = [A 



J 0.0,0-0 



free F 



(4.6) 



where N = 2L must be even for both boundary conditions. 

The fiip operator F satisfying dimi*" = dim D(u), F"^ = I, F'^ = F is defined by 



a,b 



N 

n ^a.M- 
j=0 



(4.7) 



Proceeding as for the Klein Bottle, the number of negative eigenvalues of F for the 
Mobius strip is 



2Sa 



A 



N 



k=l 



a,k 



where the edge spin is a = 1 for fixed + boundaries and a = 2 for free F boundaries. 
The matrices D{u) commute with F. In fact, using ( |2.7| ) and ( [4.5| ) we find 



FD{u)F = D'^{\ -u) = D'^{u) = D{u). 
The partition function is thus 



^Mobius/ \ 



Tr 



FDiu 



,M/2 



U 



,Af/2 



(4.9) 



(4.10) 



where the sum is over all eigenvalues. 

More generally, for other lattice models with integrable boundary conditions that 
are non-diagonal, the condition = b^ = aj^j = b^ is not satisfied and the last 
equality in ( [4.9|) does not hold. In such cases it can be shown ||^ that D{u) is similar 
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to a symmetric matrix D{u) that does commute with F. Let us define the diagonal 
matrices 

/ r. \ 1/4 AT / \ 1/4 AT 

^.,5= n^a.v -^^,1= 1^ n^.... (4-11) 

\ '-'ao / j=0 V'-'ajv / j=0 

SO that 

FAF = A-\ AF = FA~\ FA = A-^F. (4.12) 
Let us also define 

D{u) = AD{u)A'\ (4.13) 

Then, by using ( p. 71) we find 

£)^(u) = D{u) (4.14) 

so that -D('u) is real symmetric for m G R. It then also follows using the valid part of 
(^) , (|]T|) and (|]14D that F and D{u) commute 



Fb{u)F = FAD{u)A'^F = A'^FD{u)FA (4.15) 
= A-^D^{u)A = {AD{u)A-Y = if{u) = b{u) 



In this case F>(u) should be replaced by F){u) in the partition function ( [4. 101) . 



4.1 Finitized Ising partition function on the Mobius strip 



The transfer matrices D{u) for both fixed and free boundary conditions form |^ a 
commuting family D{u)D{v) = D{v)D{u) whose eigenvalues D{u) satisfy |T6[ the 
functional equation 

A cos2(^+i) 2u - sin2(^+i) 2u 

D(u)D(u + X) = (4.16) 

^ ^ ^ ^ cos4n ^ ' 

where A = n/A is the crossing parameter. For free boundary conditions the solution 
of this functional equation yields 

D{u) = 2-^ n (cosec (^|^^) + f^ksmiu^ (4.17) 

where N = 2L and yUjt = ±1 so that there are 2^ eigenvalues. The same solution 
applies for fixed + boundary conditions but with the constraint 

L 

n /ifc = 1. (4.18) 

k=l 

In this case there are 2^~^ eigenvalues in agreement with ( [4.6| ). The eigenvalues D{u) 
are thus specified by the set {/ifc}- 
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The asymptotic behaviour of the finite-size partition function Z^m of the critical 
Ising model on the Mobius strip in the limit of large N and M with the aspect ratio 
M/N fixed is given by 

"^(«) = Tt[FD{u)''/'] ~ exp [-iVM/Mik(w) - Mf\,^,{u)] Z^°^-^(g) (4.19) 

where D{u) is the double row transfer matrix, /buik(M) is the bulk free energy, /bdy(w) is 
the boundary free energy and Z^°^^^{q) is the universal conformal partition function. 
The leading corrections to each transfer matrix eigenvalue D{u) are of the form 

^ logD{u) = -iV/b,ik(w) - /bdy(w) + ^ - A - fc) sin4M + o (4.20) 

with A = 0, 1/2, 1/16 and k a non-negative integer. On the Mobius strip we define 

q = exp(- TT ^ sin4u) = IgTomsT^^ (4.21) 



where gTorus is the modular paramter on the torus ( |2.15|) . The finitized conformal 
partition function is thus 

Z^''^--mq) = j:FJ^) (4.22) 

where tt, = 0, 1, . . . labels the eigenvalues Dn of D{u) and Fn = ±1 of -F for columns 
and Dq is the largest eigenvalue with all /ifc = 1. 

Let m be the number of /ifc such that /i^ = —1. Assuming that m remains finite 

as ^ oo, we see that 

lim TV log (^^] = -TT sin 4u Y (2k - 1)S.^ _i (4.23) 



Hence 



j-j^{fc-i/2)v„-i (4.24) 



where q is given by ( [4.211 ). So the finitized partition function is 

L 



^Mobius 



iL;q) 



fc=i 

L 



(4.25) 



k=l 



where {/xj^ denotes the sequences restricted by 11^=1 /^fc = =tl- 
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From [TB|, the eigenvalues ( 4.17| ) can be written as 



D{u) = e{4^e'-y' (e"" + tan (^|^)) (4.26) 

with e = 11^=1 A^fc- As in Section p3| , these eigenvalues are determined by their zeros 
in the complex w-plane as given by ( p.62|) . Since the zeros occur in complex conjugate 
pairs we will only look at the zeros of D{u) in the upper-half plane given by 

71 71 I TT 

Uk = -yUfc - - - -logtL+i_fe mod-, A; = 1,2,...,L (4.27) 

where /2fc = Hi+i-k and tk = tan ( 4(2^+1) )• '^^^ zeros occur on certain fixed lines and 
are classified into 1-strings and 2-strings. The relevant analyticity strip is — vr/S < 
Re{u) < 37r/8 and a 1-string is a single zero Uk such that 

Re{uk) = ^ (4.28) 

o 

A 2-string is a pair of zeros {uk, u'^) with equal imaginary part and 

(Re(ufc), Re(M',)) = f-^ , ^) (4.29) 



A distribution of zeros for a typical eigenvalue D{u) for L = 8 is depicted in Figure 0. 



"2T" 



Figure 3: Zeros of a typical eigenvalue of D{u) in the period strip —A < Re(M) < 3A 
in the complex w-plane. The relevant analyticity strip is —A/2 < Re(M) < 3A/2. 



The eigenvalues are completely classified by the string contents and the relative 
orderings of the strings in the analyticity strip. For each eigenvalue D{u), there are 
L = N/2 distinct strings on the upper-half complex w-plane. Let m be the number 
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of 1-strings in the analyticity strip and n the number of 2-strings. Then the string 



contents m and n satisfy the simple (m, n) system p2|, ^ 



m + n = L. (4.30) 

Note that j = 1 labels the 1- and 2-strings closest to the real axis and for each 
eigenvalue the ordering of the imaginary parts of the strings is 

lm{uL) > Im(ML_i) > . . . > Im(-U2) > Im(Mi) > 0. (4.31) 

From ( [4.27|) , a string at lm{uk) is a 1-string if fik = —1 and a 2-string if /i^ = 1. 
Suppose that 

/ifc = -1, k = ki,k2, . . . ,km, ki < k2 < . . . < km- (4.32) 
Then we define the quantum numbers 
^ 1 

Ij = ^ -(1 + /ifc) = {# of 2-strings above given 1-string kj} (4.33) 

Alternatively, these can be defined recursively by 

Ij - + kj - kj^i = 1, J = 1, 2, . . . , m (4.34) 

with ko = 0, Iq = L — m and solution 

Ij = L — m + j — kj, j = 0, 1, . . . , m. (4.35) 

The quantum numbers / = (Ji, I2, . . . , Im) satisfy 

L-m = n>Ii>l2>...>Im>0 (4.36) 

and uniquely label all of the eigenvalues with given string content m. 

The commuting matrices D{u) and F have a common set of eigenvectors. If the 
eigenvector x associated with the eigenvalue D{u) satisfies Fx = x we say that the 
eigenvalue D{u) has parity F = +1 under the flip. Otherwise, if the eigenvector x 
satisfies Fx = —x, we say that the eigenvalue D{u) has parity F = —1 under the 
flip. For both fixed and free boundary conditions, we find the following results from 
numerics: 

1. An eigenvalue D{u) with quantum numbers (Ji, I2, . . . , Im) = (0, 0, . . . , 0), such 
that the 1-strings are all further from the real axis than the 2-strings, has parity 
F = -1-1 under the flip F. 



27 



2. An eigenvalue D{u) with quantum numbers (/i, I2, ■ ■ ■ , Im) has parity 

F = (_i)^i+^2+...+/™ (^4 37) 

under the flip F. This implies that interchanging the position of a 1-string with 
the position of an adjacent 2-string changes the sign of the parity. 

These parity properties can not be obtained by studying the eigenvalues alone. Rather 
these parity properties reflect deep properties of the eigenvectors which we do not 
study here. Nevertheless, these observations allow the flip eigenvalue F = ±1 to be 
read off directly from the eigenvalue D{u) just by looking at the locations of its zeros 
in the complex w-plane. 

For fixed + boundaries on the Mobius strip, it now follows from ( [4.25|) and our 
numerical observations that 

1^1+ k=l meven / 

L 



El 
q 



2 

m even 



m 



= Xo{L; -q) (4.38) 

[-1) 



where we have used ( ^.321 ) and the facts that m is even, /i^ = and 



E - ^) = E + 1 - - ^) = ^ + E (4-39) 

k=i ^ j=i j=i 

Similarly, for free boundaries on the Mobius strip, the finitized partition function is 

L L 



z^{L-q) = E ^ n + E ^ n ^ 



{fc-l/2)5^^,_i 



E g^E(-^)^-^'^ + ^'/' E g"^E(-?)^-^'^ 



m odd 



= Xo(L; -q) + q'^'Xy^iL; -q) (4.40) 

If we let L ^ 00 we obtain the conformal partition functions of the Ising model on the 
Mobius strip in terms of Virasoro characters. Notice also that removing the flip F in 
(f4.38| ) and ( |4.40| ) gives the flnitized cylinder partition functions Z^^(g) and Z^^^{q) 
respectively. For the Mobius strip with free boundary conditions, the even sector of 
the spectra with i? = 1 is given by Xq and the odd sector with R = —1 is given by 

-^1/2- 

Lastly we note that, in agreement with ( [4.8| ), the number of negative eigenvalues 
of F is 



lim - 

9-1 2 




Z'^\\L-q)-Z\L-q)\ = \ ' , ' (4.41) 
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where we have used 



Iff;), L odd, m even 



if PI - PI = <^ f442) 



2 



U) (lIj) 



otherwise. 



5 Summary and Discussion 

In this paper we have obtained the finitized conformai partitions of the Ising model 
on the Klein bottle and Mobius strip using Yang-Baxter techniques and the solution 
of functional equations. Changing the topology from the torus to the Klein bottle or 
the cylinder to the Mobius strip is achieved by inserting a flip operator before gluing 
the boundaries with the trace. Under the action of the flip operator, a parity F = ±1 
is assigned to the eigenvalues of the commuting transfer matrices. The task is thus 
reduced to finding the parities of the eigenvalues of the relevant transfer matrices. 
Although in this paper we have only considered the Ising model, our methods and 
observations apply more generally to solvable lattice models as we now explain. 

For the Klein bottle, it was shown analytically that all the complex eigenvalues 
of the single row transfer matrix cancel with their complex conjugates. The partition 
function for Ising model on the Klein bottle is thus given by a projection from the 
partition function on the torus 

= \Xo{q)? + \Xi/2{q)? + \Xi/M? 

^ Xo(g')+Xi/2(g')+Xi/i6(g') = Z'^'^-^q) (5.1) 

where only the left- right chiral symmetric part is retained. In fact, from numerics, 
this projection is valid for any model with L > 3 so that we observe the projection 

ziri^) = E ixA(g)r ExA(kr) = zf--{q) (5.2) 

A A 

where the sum is over the Kac table. Similarly, we have checked that the same 
mechanism works for the D4 or Potts model 

ZlT{<D = \Xo{q) + X3iq)\' + IX2/5(g) + Xv^q)]' + 2|xi/i5(g)P + 2\x2/3iq)\' 

^ xo(g') + xsiq^) + X2/5(g') + X7/5(g') + 2x1/15(5') + 2x2/3(5') (5.3) 

= z^'rii) 

For the Mobius strip, the modular parameter effectively changes from q to —q 
under the action of the flip operator which transforms the topology from the cylinder 
to the Mobius strip. More specifically, since q = exp(27rir), this is equivalent to 
saying that r i— > t + 1/2. Again we find numerically that this projection mechanism 
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applies to other models such as, for example, A4 with the (r, s) = (2, 2) boundary 
condition 

^Cyl, (2,2) (^^^j ^ XlM + XlM + XsM + xsM 

^ + xi,3(-g) + X3,i(-g) + X3,3(-g) (5.4) 

_ r^Mobius. (2,2) / \ 
— ^Ai 

where the tilde on the Virasoro characters means that the leading fractional power of 
q is replaced by the same fractional power of 

Our results for the Ising model are in agreement with 0, ^ . To compare the results 
with IfZll, we use the identities 

= g-V48 TT (1 + ^n+1/2)^ ^ ^-1/48 TT _ ^n+1/2) (5 5) 

V n=0 V V n=0 

and r = 'i^ + l/2to re-express ^fixcd{±) (t) and Zf^^^{T) in (11,12) of [0 from theta 

(2) 

and eta functions to q-product representation (Note that Zfj.^^{T) corresponds to anti- 
periodic boundary which we don't consider here). 

In 0, again we first re-express the theta functions into q-product represent aion. 
For Klein bottle, we use the fact that n^=i(l - g^"'^)(l + g") = 1 + 0(g^+^), and 
substituting the aspect ratio ^ into q, we obtain the same exact formula (|5.1| ) in 
terms of q-products ( p.20 , 2.21 , 2.21|) . For Mobius strip, we confirm the result by series 



expansion about g = up to 100 terms using Mathematica. 

Finally, we remark that it would be of interest to generalize our results on the Klein 
bottle and Mobius strip to the more general twisted boundary conditions discussed 
recently in the context of conformal field theory by Petkova and Zuber 0. 

Appendix A: Some Linear Algebra 

We list some useful properties of T{u) and F: 
1. If T{u)x = T{u)x and T{u) is real then 



T{u)x = T{u)x = T{u)x (A.l) 



so that if X is real then T{u) = T(u) and T(u) is real. 

2. If T{u) is a real normal matrix with T'^{u) = T(A — u) then the eigenvalue of 
T'^{u) = T(A — u) is T{u). This follows since, if a; is a common eigenvector of 
T{u) and T(A - u) so that T{u)x = T{u)x and T(A - u)x = T(A - u)x, then 

T(A - u)\x\'^ = x^T{\ - u)x = x^T^{u)x = {T{u)xfx = T(^|a;|2 (A.2) 



and so T(A — u) = T{u). 
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3. Suppose further that Tj and Tj are eigenvalues of T = T{uq) and Tj 7^ T,-, then 
the corresponding eigenvectors Xi and orthogonal. This follows since 



(A.3) 



which implies (Tj — Tj)x\xj = and xjajj = 0. 
4. For any eigenvector a; of T{u), Fx is also an eigenvector of T{u) with the 



corresponding eigenvalue T(u): 



T{u) (Fx) = {T{u)F)x = {FT{X - u)) x = T{u)Fx 



(A.4) 



5. By ( p.l6| ) and the normality of T{u), F and D(u) = T(A — u)T{u) commute: 
F(T(A - u)T{u)) = T{u)FT{u) = T{u)T{\ - u)F = (T(A - u)T{u))F (A.5) 

Appendix B: Euler Maclaurin 

For the single row transfer matrix T{u), we compute the scaling limit of the largest 
eigenvalues in the R = +1 and R = —1 sectors by applying the midpoint and endpoint 
Euler-Maclaurin formulas 

E /(a + {2k - \)h) = -J^ fit) dt - - /'(a)] + Qih^) (B.l) 

m 1 1 h 

where 6 = a + mh. These formulas are valid if f{t) is twice differentiable on the closed 
interval [a,b]. 

In the R = +1 sector, the logarithm of the largest eigenvalue can be written as 



1 1 ^ 

logTo(M) = -(1 — L) log2 H — ^ log sin 4m + esc 



T(2fc-1) 

2L 



k=l 

Suppressing the dependence on u, we define the function 

/(t) = log t{n - t) (sin 4m + csct) 

so that 



logTo(M) = -(l-L)log2 + ^/ 
^ k=i 



7r(2fc-l) 
2L 



fc=i 



7r(2fc-l) 
2L 



(B.3) 



(B.4) 



(B.5) 
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We define / in this way to remove the singularities of csct at t = 0, vr. Applying 
the midpoint Euler-Maclaurin formula and the asymptotic expansion of the gamma 
function 



log T{z) = - I) log ^ - z + I log 27r + — + 0{z-^) (B.6) 



1 

we obtain 

N 1 ^ . ^ , 71 , /I 



logTo(M) = — log-(sin4M + csctjc/t + ^sin4M + o(^— j (B.7) 



which agrees with (|2.16|) with c=l/2, A = A = /c = /c = and bulk free energy 

1 r. 1 



/bulk = ~ ^ 2 ^ ^ ^^'^"^ 



Similarly, using the endpoint Euler-Maclaurin formula in the R = —1 sector, the 
logarithm of the largest eigenvalue is 

logTi(u) = - logL + — — log2 + - ^ log [sin 4m + esc 

fe=i 



— / log - 1 sin 4m + CSC t (it sin 4m + — (B.9) 

47rio 2V / 6A^ \nJ ^ ' 



with /bulk defined by ( [B.8| ) and A = 1/16. Hence we conclude that 

log [Ti(m)/To(m)] = sin 4m + o (1) (B.IO) 



and we obtain the factor in (|3.87|) 
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